ON STABLE SELF SIMILAR BLOW UP FOR EQUI VARIANT WAVE MAPS: 

THE LINEARIZED PROBLEM 



ROLAND DONNINGER, BIRGIT SCHORKHUBER, AND PETER C. AICHELBURG 



Abstract. We consider co-rotational wave maps from (3 + 1) Minkowski space into the three- 
sphere. This is an energy supercritical model which is known to exhibit finite time blow up via 
self-similar solutions. The ground state self-similar solution /o is known in closed form and based 
on numerics, it is supposed to describe the generic blow up behavior of the system. In this paper we 
develop a rigorous linear perturbation theory around /o . This is an indispensable prerequisite for the 
study of nonlinear stability of the self-similar blow up which is conducted in the companion paper 
In particular, we prove that /o is linearly stable if it is mode stable. Furthermore, concerning 
the mode stability problem, we prove new results that exclude the existence of unstable eigenvalues 
with large imaginary parts and also, with real parts larger than i. The remaining compact region 
is well-studied numerically and all available results strongly suggest the nonexistence of unstable 
modes. 



1. Introduction 

Let {M,g) be a Lorentzian manifold and {N,h) a Riemannian manifold with metrics g and h, 
respectively. Consider the geometric action functional 

S{U) = [ tlg{u*h), 

Jai 

for a function u : M ^ N, i.e., S{u) is the integral over the trace (with respect to g) of the pullback 
u*h. In abstract index notation, this expression reads 

Siu) = [ g''''id^u'')id,u')habOU 
Jm 

and formally, in local coordinates (x'^), the Euler-Lagrange equations associated to this action are 
given by 

(1.1) Ugu'^ix) + g^''{x)Tt,{u{x))d^u\x)d^.u'{x) = 

where Einstein's summation convention is in force. The indices /i, v and a, 6, c take the values 
0, 1, ... , dim M — 1 and 1,2,... , dim A^, respectively, and denotes the Laplace-Beltrami operator 
on M. Furthermore, F^^ are the Christoffel symbols associated to the metric h on the target 
manifold A^. Eq. (II. ip is known as the wave maps equation in the intrinsic formulation. In order to 
develop some intuition, let us consider simple examples. If = M" with the standard Euclidean 



The first author is an Erwin Schrodinger Fellow of the FWF (Austrian Science Fund) Project No. J2843. He wants 
to thank Wilhelm Schlag for many interesting discussions and the ESI Vienna for hospitality during the workshop 
"Quantitative studies of nonlinear wave phenomena". 

The second author acknowledges partial support from the program "Innovative Ideas" of Vienna University of 
Technology; the Austria- Croatia Project HR 10/2010; the Austria-France Project FR 07/2010; and the Austria- 
Spain Project ES 08/2010 of the Austrian Exchange Service (OAD). 

The third author acknowledges partial support from the FUNDACION FEDERICO and the Erwin Schrodinger 
Institute for Mathematical Physics, Vienna, for hosting a program on "Quantitative studies of nonlinear wave 
phenomena" . 

1 



metric and M is (m + l)-dimensional Minkowski space with coordinates (t, ), then 

Eq. (jl.ip reduces to 

Ufj - A]KmU" = 0, 

the ordinary Hnear wave equation for a vector-valued function. As a consequence, if the target is 
curved, the wave maps equation constitutes a nonhnear generahzation of the wave equation and its 
geometric origin makes this generahzation very natural. Indeed, consider for instance M = (0,1), 
an open interval, and N any Riemannian manifold. Then the corresponding wave maps equation 
reads 

^W + rw«(*))^(t)^W = o 

which is nothing but the geodesic equation on N. These simple examples already indicate that the 
wave maps functional is a rich source for interesting nonlinear equations. 

In physics, wave maps first appeared in the 1960s as the SU(2) sigma model [17J which corre- 
sponds to mappings from physical (that is, (3 + l)-dimensional) Minkowski space to the three- 
sphere. In fact, the present work deals exactly with this model subject to a certain symmetry 
reduction, see below. Let us mention at least two more applications of wave maps in physics. Due 
to its geometric nature, the wave maps equation is frequently used as a toy model for the much 
more involved Einstein equations of general relativity, in particular in the context of critical gravita- 
tional collapse. Furthermore, self-gravitating wave maps are used in general relativity to study the 
behavior of nonlinear matter fields coupled to Einstein's equations. These works include the study 
of singularity formation ( "hairy black holes" , "naked singularities" ) , solitons with nonzero cosmo- 
logical constant and, more recently, long-time asymptotics ("tails") of global solutions (see e.g., 
[9], [28], [3] and references therein). We also refer the reader to (31] for more (possible) applications 
of wave maps in physics. 

From now on we restrict ourselves to wave maps on Minkowski space, i.e., we assume that 
M = W^~^^ with the standard Minkowski metric {qhu) = 'iiag(— 1, 1, . . . , 1) and coordinates 
{t, x^,x'^, . . . , x"'). In this case, Eq. (|l.ip reads 

(1.2) <t(t, x) - Avnu^it, x) = g>"'TtMt, x))d^u\t, x)d^u%t, x) 

and it is natural to consider the Cauchy problem, that is, one prescribes initial data (m'^(0, •), (0, •)) 
and studies the future development. The Cauchy problem for the wave maps equation has attracted 
a lot of interest in the past 20 years. Eq. ()1.2p is a system of semilinear wave equations and the 
standard theory immediately yields local well-posedness in the Sobolev space H^(W^) x H^~^{W^) 
provided that s > y -|- 1, see e.g., [22]. However, by exploiting the null structure of the wave maps 
nonlinearity, it is possible to improve this result to s > ^ ( "Klainerman-Machedon theory"), see 
the survey [22] and references therein. The wave maps equation ()1.2p is invariant under the scaling 
u(t, x) I— )• ux{t, x) := u{Xt, Ax) for a A > and there exists a conserved energy given by 



dx 



hab{u{t, x))dtu"-{t, x)dtu''{t, x) + hab{u{t, x)) ^ diu"-{t, x)diu^{t, x) 

1=1 

which scales like E{ux) = A^~™£'(ti). This shows that the wave maps equation is energy subcritical 
if m = 1, critical if m = 2 and supercritical if m > 3. Based on a well-known heuristic principle 
one expects large data global well-posedness in the energy subcritical case and finite time blow up 
for energy supercritical equations. Indeed, large data global well-posedness in the subcritical case 
has been proved in [20]. Also, small data global well-posedness in dimensions m > 2 for a large 
class of targets has been estabhshed in the last ten years, see e.g., [50], [51], [36], [17], pT^, [38], 
[52] , [26] , [27] , [33] , [32] . We refer the reader to [23] for a survey on these results and a detailed list 
of references. In the energy critical dimension m = 2, the questions of large data global existence 
and blow up are related to the geometry of the target. For spherical targets (and also more general 
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surfaces of revolution), blow up solutions have been constructed in [25], [35] and, very recently, 
[5], [33]. Again, we refer the reader to p3] for more references and earlier results related to this. 
On the other hand, large data global well-posedness results for the energy critical case have been 
obtained in [Mj as well as [38], [39] and in the series [35] . 

In the present paper we study the energy supercritical case. More precisely, we consider wave 
maps from physical Minkowski space to the three-sphere, the original SU(2) sigma model of particle 
physics. Furthermore, we impose a restrictive symmetry assumption: we require the wave map to 
be co-rotational. In order to explain what this means, choose standard spherical coordinates 
{t,r,9,ip) on Minkowski space and hyperspherical coordinates {tpjQ,^) on the three-sphere, i.e., 
the respective metrics are given by 

g = -dt^ + dr^ + r^idO'^ + sin^ Odip'^) 

and 

h = di;^ + sin^ ^{de'^ + sin^ Qd<l>'^). 

In these coordinates, a mapping u : M ^ N is described by the three functions tp{t,r,9,{p), 
@{t, r, 0, if) and r, 9, ip). The mapping u is said to be co-rotational if Q and ^> are trivial in the 
sense that Q{t, r, 9, ip) = 9, r, 9,ip) = ip and ip is independent of 0, (p. The wave maps equation 
for co-rotational maps reduces to the single semilinear wave equation 

/ , , 2 , sin(2V') 

(1.3) ^tt - Ar - -A + — ^ = 

and requiring regularity (e.g., C^) of ijj at the origin yields the asymptotic behavior il){t,r) = 0(r) 
for any t as r — )• 0+ and in particular we have the boundary condition -(/^(t, 0) = for all t. However, 
note carefully that smoothness of at r = does not imply that iprit^Q) has to vanish (observe 
the special cancellation by Taylor expansion). 

A considerable amount of the wave maps literature is devoted to the study of the Cauchy problem 
in the presence of symmetries such as equivariance or spherical symmetry, see e.g., |7], [8], [33], [32], 
[33], [E], [31]) [37]. For the model under investigation, the papers [30] and [39] are of particular 
interest. In [40j . small data global well-posedness in a sufficiently high Sobolev space is established 
whereas in [39], many different aspects of the Cauchy problem are considered, in particular the 
minimal regularity requirements for local well-posedness. It is well-known that Eq. (jl.3p exhibits 
self-similar finite time blow up. A solution ijF of (II. 3p is said to be self-similar if it is of the 
form ^'^{t,r) = /{jT—j:) for a constant T > and a smooth function / 0. The existence of a 
smooth self-similar solution /o for Eq. ()1.3p has been proved in [36] by variational techniques and 
independently, it has been found in closed form [53], see also [6]. We refer to /o as the fundamental 
or ground state self-similar solution and reserve the symbol ip'^ for this solution, i.e.. 



2 arctan \ rp ^ 



By exploiting finite speed of propagation, /o can be used to construct a perfectly smooth solution 
of Eq. (jl.3p with compactly supported initial data that breaks down at t = T. However, a natural 
question to ask is how generic this break down is. Does it only happen for very special initial 
data or is there an "open" set of data that lead to this type of self-similar blow up? Numerical 
experiments [3] indicate that the latter is true, i.e., the blow up described by /o is conjectured to 
be stable. More precisely, one observes that the future development of sufficiently large generic 
initial data converges to ilP- in the backward lightcone of the spacetime point (T, 0) as t — t- T— . 

It is worth mentioning that the blow up behavior of the corresponding energy critical model of 
equivariant wave maps on (2 + 1) Minkowski space is fundamentally different. Due to a result of 



"'^A priori, a self-similar solution is of the form Tp{t,r) = /(f) but one may immediately apply the time translation 
and reflection symmetries of the equation to obtain the one-parameter family ip'^{t,r) — /(j^). 
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Struwe [l3], it is known that the blow up in the energy critical case cannot be self-similar. In this 
respect it is interesting to note that Struwe's result does not imply that blow up actually happens. 
However, as already mentioned, blow up solutions for equi variant wave maps from (2+1) Minkowski 
space to the two-sphere have been constructed in [25], [35], [3^ . 

1.1. Outline of the paper. In the present paper we develop a functional analytic linear perturba- 
tion theory around the self-similar solution ip'^ . Since ip'^ is time-dependent, we introduce adapted 
coordinates (r, p) ("similarity variables") such that ijF becomes independent of the new time co- 
ordinate r. The coordinates (r, p) cover the backward lightcone of the blow up point (f, r) = (T, 0) 
which is mapped to r = oo. The linearization of Eq. (II. Sp around the solution ifp" in the adapted 
coordinates is written as a system of the form 



where L is a linear spatial differential operator which is realized as an unbounded operator on a 
suitable Hilbert space. However, the linearized operator L is not self-adjoint (in fact, not even nor- 
mal) and thus, the analysis of this equation is highly nontrivial since one cannot resort to standard 
methods from self-adjoint spectral theory. Consequently, we apply semigroup theory to study the 
linearized time evolution of perturbations of ijF . In Sec. [21 we present a heuristic discussion that 
motivates the choice of the underlying Hilbert space. Then, in Sec. O we review known theoretical 
and numerical results on the point spectrum of L, that is, we consider mode solutions ^(t) = e'^'^u 
of Eq. (jl.4p and prove a new result that excludes unstable eigenvalues with real parts larger than 
^. In Sec. m we show that L generates a strongly continuous one-parameter semigroup S{t). This 
yields the well-posedness of the Cauchy problem for Eq. (jl.4p and provides the rigorous spectral 
theoretic basis for the mode analysis in Sec. [3l It turns out that the time translation symmetry of 
the wave maps equation induces a single unstable eigenvalue in the spectrum of L. This introduces 
an additional difficulty since the instability is "artificial" and one is actually interested in the flow 
"modulo this instability" . In order to make this precise, we construct an appropriate spectral pro- 
jection that decomposes the underlying Hilbert space in a stable and an unstable part. We show 
that the unstable subspace is spanned by the single eigenvector associated to the aforementioned 
unstable eigenvalue of the linearized operator L. Finally, we prove an appropriate growth bound 
for the linear time evolution S{t) restricted to the stable subspace which shows that iji^ is indeed 
linearly stable, provided that there are no unstable modes. As a consequence, we obtain a thorough 
understanding of the linear stability problem for ■0-^. Remarkably, we identify many features that 
are known from analogous self-adjoint problems. As a by-product, we also show that unstable 
eigenvalues with large imaginary parts do not occur. This result justifies some of the numerical 
methods that have been employed to study mode stability of ■0"^, in particular the shooting method. 
The results of this paper are of fundamental importance for the study of the full nonlinear stability 
of tlP" which is conducted in jllj . 

1.2. Notations and conventions. For Banach spaces X^Y we denote by B{X,Y) the Banach 
space of bounded linear operators from X to Y . As usual, we write B{X) \i X = Y . Vectors 
are denoted by bold letters and the individual components are numbered by lower indices, e.g., 
u = (tii,ii2). We do not distinguish between row and column vectors. Furthermore, we denote the 
spectrum and resolvent set of a linear operator A by (j{A) and p{A), respectively. In particular, 
we write (y,p[A) for the point spectrum, i.e., the set of all eigenvalues. For A E p[A) we set 
-Ra(A) := (A — A)~^, i.e., Ra is the resolvent of A. Finally, for a, 6 G M we use the notation a < 6 if 
there exists a constant c > such that a < cb and we write a ~ 6 if a < 6 and b < a. The symbol 
~ is reserved for asymptotic equality. Also, the letter C (possibly with indices) denotes a generic 
nonnegative constant which is not supposed to have the same value at each occurrence. 



(1.4) 
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2. Self-similar blow up 



Our strategy is to treat the nonlinearity in Eq. ()1.3p as a perturbation of a suitable "free" 
equation. Note that the nonhnearity is singular at r = and thus, we have to split off the singular 
part first and assign it to the differential operator, i.e., the free equation. In order to do so, we 
simply rewrite Eq. ()1.3p as 

. X , , 2 , 2 sin(2V') - 211) 

Formally, the nonlinear term is now regular at r = since we assume il>{t,r) = 0{r) as 

r —7- 0+. Our goal is to study Eq. (|2.5|) in the backward lightcone 

CT:={{t,r):tG{0,T),rG[0,T-t]} 

of the blow up point {t,r) = (T,0). The conserved energy of the free equation 

2 2 

(2.6) Iptt - Iprr - Z'^r + — = 

is given by 

1 
2 







as follows immediately by multiplying the equation by r'^Tpt and integrating by parts. Naively one 
might expect this energy to yield a natural norm suitable for the study of the stability of ip'^ . We 
give an argument why this is not the case. Consider the local energy in the backward lightcone Ct 
given by 

Mi) ■= o / (rV?(t,r) +rV'(i,0 + 2V''(t,r)) dr. 
^ Jo 

By energy conservation, one has the estimate E^{t) < 1 and for a generic energy solution living on 
Ct this cannot be improved to obtain polynomial decay of the form E^{t) < (T — ty as t ^ T— 
for some constant 7 > 0. To see this, consider the function 



V^,(t,r) = (T-t)-2+-n4^ 



where 



for arbitrary e > which is a solution of the free equation inside the lightcone Ct- Its local energy 
satisfies E^^{t) ~ (T — t)^^ and, since e > is arbitrary, this shows that there cannot exist a 7 > 
such that E^it) < (T — t)'^ as t ^ T— for generic solutions -0 of the free equation. On the other 
hand, for the self-similar function ijp" we have E^T{t) < {T — t) and thus, the local energy of ilF 
decays faster than the local energy of a generic free solution. This, of course, renders a perturbative 
approach hopeless from the very beginning. As a consequence, the local energy does not yield a 
suitable norm for studying the blow up of ilP" because it does not blow up. Obviously, the same 
is true for the energy of the full nonlinear equation (II. 3p . Thus, in order to make a perturbative 
approach feasible, we need to find a norm || • || which is derived from a conserved quantity of the free 
equation and \\{'4F{t, ■),ipj'(t, •))|| — )• 00 as t — )■ T— . The idea is obvious: we have to consider the 
gradient instead of the function itself since hitting the self-similar solution ip'^ with brings out 
the singular factor {T — t)^^. Thus, we need a convenient equation for ipr- However, differentiating 
Eq. (12. 6p directly yields a true third-order equation for ipr due to the term which is not 
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desirable. Hence, we have to remove this term first by transforming to a new dependent variable 
ip of the form ip{t,r) = r^'ijj{t,r) for a suitable k. We obtain 

~ ~ 2{k + l) ~ {k-l){k + 2) ~ 

i^tt - At - -^r - = 

and this shows that there are two possibilities, either A; = 1 or A; = —2. Let us investigate those 
two possible paths. At first glance one might expect fc = 1 to be the more appealing choice since 
it leads to the radial 5-dimensional wave equation 

4 . 

- V'rr V'r = 

r 

and differentiating with respect to r yields 

4 , 4 . 

Multiplying by r'^iprt and integrating by parts we obtain the conserved quantity 
(2.7) ^ (r^Atit, rf + r^rr{t^ rf + Ar^r{t, rf) dr. 

This suggests to use ipt and V'r as evolution variables. However, there is a serious disadvantage of 
this choice: we do not have a boundary condition for ^^(t,r) (recall that V' = ■f and ijj{t,r) = 0{r) 
as r ^ 0+). As a consequence, at a fixed time slice t = const, we cannot reconstruct the orginal 
field ifj from the proposed evolution variables 'i/'t, ipr- This shows that the full equation (|1.3p cannot 
be written in terms of the variables ipt, ipr- Another way to put this is to note that the prospective 
local energy "norm" derived from the quantity ()2.7p is in fact only a seminorm on the space of 
functions we are interested in which is too weak for our purposes. Consequently, we dismiss the 
choice k = 1 and follow the second possible path. 
Choosing k = —2 we obtain the equation 

2 ~ 

Iptt - Iprr + -tpr = 0. 

r 

Note carefully that this is not a radial wave equation since the first order term has the wrong sign! 
Differentiating with respect to r yields 

2 ~ 2 ~ 

In order to obtain a conserved quantity ("energy") for this equation, one notes that ^^rr — fV'rr = 
r'^dr{r~'^iprr), multiplies by ^iprt and integrates by parts. However, the resulting expression is not 
positive definite due to the wrong sign of the first-order term. Thus, we make another transforma- 
tion tpritjr) =: r^'ijj{t,r) to obtain 

7 ? k-l ^ (k-l)(k-2) ; 

iPtt - i^rr - A - ^ = 

and by choosing k = 1 we end up with the one-dimensional wave equation 

At -1prr = 

on the half-line r > 0. In terms of the original field we have 



ijj{t, r) = rAit, r) + 2ip{t, r) 
6 



and this implies the boundary condition 1^(1,0) = for all t. As a consequence, we obtain the 
conserved quantity 

(2.8) (7Pt{t,rf +Mt,rf)dr 



and this suggests to use ipt and ^pr as evolution variables. This is indeed possible now since the 
boundary condition V'(*)0) = allows us to express ip in terms of V'r as 

ip{t,r) = / 'i/jr{t,r')dr' 







and as a consequence, for the original field ip, we obtain 



Tp{t,r) = 4 / r''>p{t,r')dr' = \ I r' I il^rit, s)dsdr' . 



J.2 I ' ^2 



Now observe further that for the self-similar solution 

ip'^(t, r) = rifl (t, r) + lijF it^ r) = ^-^ + 4 arctan p 

1 + 

we have 



r) = {T- t)-''^-^^ =■■ {T - tr'gip) 

and 

i^f{t,r) = {T-t)-'pg{p) 

where p = YZf. This implies 

^'"^'^ {4^1{t,rf + ^Pj{t,rf) dr={T- t)-^ [(t^)^ (t^)' + 9 (t^)'] dr 

= {T-t)-' (\ph{p?+9{p?)dp 
Jo 

as t — )• T— and thus, the local "energy" norm derived from the conserved quantity (12. Sp blows up. 
Therefore, it is a promising candidate for a suitable norm to study the stability of tp'^ . In this 
paper we develop the linearized perturbation theory about tp'^ in the topology given by this higher 
energy norm. 

3. Mode stability 

3.1. Similarity coordinates. In this section we study so-called mode solutions of the linearized 
equation around the fundamental self-similar solution tp'^ . We linearize Eq. ()1.3p by inserting the 
ansatz ip = ip'^ + ip and neglecting terms of order ip"^ or higher. This yields the linearized equation 

, , 2 2cos(2'0'^) 

(3.9) iptt - '^rr H 5 = m Ct 

with appropriate initial data. The concrete form of the data is irrelevant for the mode analysis and 
therefore, we omit them in this section. Since the solution is self-similar, it is natural to switch 
to similarity coordinates (r, p) defined by r := — log(r — t), p := ij^. Observe that this change of 
variables transforms the cone Ct to the infinite cylinder 



Zt:={(t,p) :T>-logT,pG [0,1]}, 
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Figure 3.1. The backward lightcone Ct in it,r) and (t, p) coordinates. We have 
also included some lines of constant p (dashed). 



see Fig. 13.11 In particular, the blow up of V'^ takes place as r goes to infinity. According to the 
chain rule, the derivatives transform as dt = e.^{dr + pdp) and dr = e^dp. In similarity coordinates 
(r, p) , the evolution equation (j3.9p reads 

ri2 



(3.10) 
where 
(3.11) 



+ 0^ + 2p(t)rp - (1 - 



1 



+ 



V{p) 







p p^ 

{t,p) := ip{T — e~'^,e^'^p) and the "potential" V is given by 

2(1 - 6/92 + /) 



V{p) = 2cos(4arctan(p)) 



2\2 



The first step in a stability analysis of such an equation is to look for mode solutions, i.e., solutions 
of the form </>(r, p; A) = e'^'^u{p, A) for a complex number A and a sufficiently regular function «(•, A). 
The nonexistence of mode solutions with ReA > is a necessary (though, in general, not sufficient) 
requirement for linear stability of the fundamental self-similar solution ^ifF . 

3.2. The eigenvalue equation. In order to decide whether there are unstable mode solutions or 
not, we insert the ansatz </>(r, p; A) = e^'^u{p, A) into Eq. (j3.10p to derive the ordinary differential 
equation 

1 - -2 



(3.12) 



{l-p')u"{p,X) 



■^u'{p, A) + 2Apn'(p, A) + ^u(/). A) + A(A + l)n(p. A) = 



p p- 

which constitutes a generalized eigenvalue problem for the function ti(-,A) where ' always means 
derivative with respect to p. Unfortunately, due to the complexity of the potential V , this equation 
cannot be solved explicitly and it seems to be a hard problem to exclude unstable mode solutions 
rigorously. Indeed, there are only partial results available, see below. The rest of this section is 
devoted to the analysis of Eq. ()3.12p . 

We start off with some basic observations. Eq. (j3.12p has two regular singular points at /? = 
and p = 1 (which is the location of the lightcone). Furthermore, a brief inspection of the potential 
V shows that additionally there are regular singular points outside the "physical" domain at p = 
— 1, ibi as well as p = 00. By the change of variables p 1— p^ one can reduce the number of regular 



singular points from six to four which shows that the general solution is given in terms of Heun's 
functions. This, however, is not of much use for us since there is no rigorous theory of Heun's 
functions and most of the results are based on some sort of numerics. Nevertheless, we can at least 
get information on the asymptotic behavior of the solutions around the singular points p = and 
p = 1 by applying Frobenius' method (see e.g., [30], §6.2). The indices at p = are {1, —2} and at 
p = 1 we have {0, 1 — A} which shows that analytic solutions behave like p around and like 1 at 
p = 1. However, note carefully that the nonanalytic solution around p = 1 behaves like (1 — p)^~'^ 
and it therefore becomes more and more regular as ReA decreases. It is not clear at this stage (and 
in fact it is false) that only analytic solutions of Eq. (j3.12p are relevant. The precise regularity 
requirements can only be derived once one has a well-posed initial value formulation of Eq. (jS.lOp 
and a well-defined spectral problem. This will be postponed to Sec. [Hand for the moment we only 
look for analytic solutions of Eq. (j3.12p . Thus, we say that A E C is an eigenvalue iff Eq. (|3.12p has 
a nontrivial analytic solution n(p. A). Note also that it is sufficient to consider smooth solutions 
(i.e., solutions in C°°[0, 1]) since any smooth solution of Eq. (|3.12p is automatically analytic by 
Frobenius' method. 

3.3. The gauge mode. By a direct computation one immediately checks that the function g{p) = 
j^j^ is a solution of Eq. (j3.12p with A = 1 and thus, there exists an unstable mode. However, it is 

well-known that this instability emerges from the fact that is a family of solutions (depending 
on the free parameter T) rather than a single one. Hence, the time translation symmetry of the 
original wave map equation ()1.3p is reflected by this unstable mode of the linearized operator and 
therefore, the mode g is sometimes called a gauge mode. As a consequence, we really want to study 
stability "modulo" this gauge instability. We will make this precise in Sec. HI For the moment we 
just ignore this gauge mode and look for other unstable modes. Note also that there are no other 
"gauge modes" since all the remaining symmetries of the equation are eliminated by our restriction 
to co-rotational wave maps (of course, there is still the scaling symmetry but, when acting on a 
self-similar solution, scaling and time translation are equivalent). 

3.4. Known results. The eigenvalue equation (j3.12p has been extensively studied by several au- 
thors. First, there are a number of numerical works where essentially three different methods have 
been applied to show that there are no unstable modes 0. Probably the easiest thing to do is 
to employ a so-called shooting and matching technique where one integrates the equation (j3.12p 
away from the singular points by using a standard numerical method like Runge-Kutta or variants 
thereof. One then tries to match the two solutions at p = ^. Such a numerical analysis has been 
performed in [10] . |12j . [3], [29]. Special care has to be taken since one has to solve singular initial 
value problems but this can be done in a satisfactory way. A more serious disadvantage of this 
method is that one has to have good initial guesses and therefore, it is in principle possible to 
"miss" eigenvalues. That is why one would like to have additional (theoretical) arguments that 
there are no eigenvalues far away from the real axis. In fact, we will prove such a statement, see 
Lemma 14.231 below. Second, one can numerically integrate the time evolution equation ()3.10p and 
remove the gauge instability by a suitable projection [lOj, ti2j. With this approach one cannot 
"miss" eigenvalues with large imaginary parts but it is numerically more delicate since one has to 
solve a PDF instead of an ODE. Finally, it is possible to rewrite the problem (I3.12P as a difference 
equation and use continued fractions expansions which eventually leads to a root finding problem 
for a transcendental function. This is a highly accurate method and it has been successfully ap- 
plied to Eq. (I3.12[) in [2j. The outcome of all these numerical studies is that the eigenvalue with 
the largest real part (apart from the gauge mode) is A ~ —0.542466 and this strongly suggests the 
mode stability of the solution ip'^ . 



'We always mean unstable modes apart from the gauge mode. 
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As far as rigorous arguments are concerned, there are only partial results. In [14j . Eq. (j3.12p is 
transformed into a self-adjoint Sturm-Liouville form and this approach can be used to show that 
there are no unstable eigenvalues A with ReA > 1. Furthermore, by ODE techniques it has been 
proved that there are no real unstable eigenvalues except for A = 1 [13]. However, there is no result 
so far that excludes unstable eigenvalues A with < ReA < 1 and ImA ^ 0. In what follows we 
push the boundary a bit further and prove a result that excludes eigenvalues A with ReA > ^. 

3.5. The dual eigenvalue problem. In [14] there is given a factorization of the linear operator 
which can be used to obtain a "dual" eigenvalue problem which gives important new insights into 
the problem of mode stability. 

Lemma 3.1. Let n(-,A) G C°°[0, 1] be a nontrivial solution of Eq. (j3.12p . Then either A = 1 and 
u{p, 1) = for a constant c G C (i.e., u{-, A) is a multiple of the gauge mode) or the equation 

(3.13) - (1 - p^)v"{p, A) - 2^^^^v'{p, A) + 2Xpv'{p, A) + ^^v{p, A) + A(A + l)v{p, A) = 0, 

p p^ 

with the potential 

1 + /O^ 

has a nontrivial solution v{-,X) £ C°°[0, 1]. 

Proof. Let u{-,X) E C°°[0, 1] be a nontrivial solution of Eq. (j3.12p . We define a new variable 
u{p,X) := (1 - />2)V2^(p^ A). Then n(.,A) satisfies 



-(1 - pru"{p, A) - 2iL^n'(/,, A) + P\ ip, A) = -(1 - X)Mp, A) 

p p^ 



and this equation can be written as 

(3.14) /3/3n(-,A) = -(l-A)2^2(.,A) 

where the (formal) differential operators f3 and f3 are given by 

dp p{l + p^) 

and 

d 3-/32 



13 :=-{!- p^)- 



'dp p{l + p^y 

see [11]. The solutions of the first-order ODE (3u = can be given explicitly and the one- 
dimensional solution space is spanned by the function p 1— )• ^Y+p^^ • Thus, if /3n(-,A) = then 
Eq. (|3.14p implies A = 1 and n(-,A) is a multiple of the gauge mode. Suppose f3u{-,X) ^ 0. 
Applying the operator f3 to Eq. (|3.14p yields 

(3.15) /3/3/3n(.,A) = -(l-A)2/3u(.,A) 

and this is a second order ODE for the nonzero function v{-,X) := /3n(-,A). The operator /3 is 
singular at p = but due to the asymptotics u{p, A) ~ p as p — )• 0+ (Frobenius' method) we 
still have v{-,X) G C°°[0, 1). Furthermore, as p — ?■ 1—, ^(p. A) behaves like (1 — p)^^"^- Explicitly, 
Eq. (f3l3D reads 



-(1 - prv"ip, A) - 2^^-^v'ip, A) + (1 P"^\(P^ P% ^p, A) = -(1 - X)Mp, A) 



p p^ 



with the new potential V{p) = ■ Setting v{p,X) := (1 — p^) ^^'^v{p,X) we obtain v{-,X) S 

C~[0, 1] and v{-, A) satisfies Eq. (IXTSll . □ 
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3.6. Nonexistence of eigenvalues with ReA > ^. Based on Lemma |3. II we can now prove that 
the gauge mode is the only eigenmode with ReA > |. 

Proposition 3.2. Let u{-,X) G C°°[0, 1] be a nontrivial solution of Eq. (|3.12p with ReA > |. Then 
A = 1 and there exists a (nonzero) constant c G C such that u{p, 1) = j^Tpi- 

Proof. Let u{-,X) G C°°[0, 1] be a nonzero solution of Eq. (j3.12p with Re > ^ and assume that 
A 7^ 1. According to Lemma |3 . 1 1 ther e exists a nonzero v{-, A) G C°°[0, 1] such that 



(1 - p^)v"{p, A) - 2 —v'{p, A) + 2\pv'{p, A) + -^v{p, A) + A(A + l)v{p, A) = 0. 



P 



p- 



Thus, /9;A) := e v{p,X) is a smooth solution of the evolution equation 

Xtt + Xt + '^PXrp - (1 - P^)Xpp - '2- —Xp + -^-^X = 

p P'^ 

and transforming back to the physical coordinates (t, r) this simply reads 

Xtt - Xrr Xr H 5 X = 0. 

Motivated by the energy of this equation, consider the positive smooth function 



1 



T-t 



2 JO 

Differentiation yields 



r'\Xt{t,r;X)\'+r'\Xr{t,r;X)\' + V 



T-t 



\x{t,r;X)\' 



dr. 



F'Jt) = Re 



T-t 



r^XttXt + r XrXrt + V 



T-t 



XXt + :,V'{^^j^^^^\x\'dr 



<o 



- ^(^ - t)'\xtit, T - t; A)|2 - l(r - T - A)|2 - ^y(l)|x(t, T-t; X)\^ 

.T-t 

<Re [drir'^Xr)xt + t'^XtXh] dr 

Jo 

-^iT- tf\xtit,T- t; A)|2 - i(r - tflxrit, T - t; A)|2 - ^V{l)\x{t, T-t; X)\^ 
<-|x(t,T-t;A)p 

by integration by parts and Cauchy's inequality. Inserting x{'t,T — t; X) = {T — and 
integrating from to t for t < T yields 



F^{t) < -\v{l)\' [\t - tr^'^'^'dt' + F^(0). 
Jo 



However, from the asymptotic behavior of v we know that v{l) ^ and therefore we have a 
contradiction since Fy. is positive but the right-hand side of the above inequality goes to — oo as 
t —7- T— (recall that we assume ReA > □ 
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3.7. Reformulation as a resonance problem. We conclude this section by showing that the 
mode stabihty problem is in fact equivalent to finding resonances of the 1 = 2 Schrodinger operator 
with a potential that has inverse square decay towards infinity. We formulate the problem for the 
"dual" eigenvalue equation 



-(1 - p^)v"{p, A) - 2^—^v'{p, A) + 2\pv'{p, A) + ^^v{p, A) + A(A + l)v{p, A) = 



p p^ 



with potential V{p) = \^pi ■ We are interested in smooth solutions and according to numerics 
we expect the eigenvalue with the largest real part to be A ~ —0.54. Now suppose f(-,A) is an 
eigenfunction, i.e., a smooth solution to the above equation. By Frobenius' method we have the 
asymptotic behavior v{p, X) ^ p^ as p ^ 0+ and v{p, X) ~ 1 as p — )• 1—. Defining a new variable v 
by 

v{x, p) := (1 — tanh^ x)'^^'^ sinhx ti(tanhx, 1 — ip), 
the eigenvalue problem transforms into 

(3.16) - v"(x, A) + -^v(x, A) + Q(x)v(x, A) = p'^v(x, A) 

on the half-line x > with the potential 

y(tanh2;) 6 
smh X x-^ 

The outgoing Jost solution fj^{-,p) is defined as the solution of Eq. (|3.16|) that behaves like 
f+{x,p) ~ e*'^^' as X — )■ oo. Furthermore, we denote by fo{-,p) the solution of Eq. ()3.16p which is 
square-integrable near x = 0. The eigenfunction v{-,p) has the asymptotic behavior v{x,p) ~ x^ 
as X — )• 0+ and v{x,p) ~ e^^^ as x — )• oo. Thus, finding eigenvalues amounts to matching the 
outgoing Jost solution f^{-,p) to fo{-,p) or, in other words, A is an eigenvalue iff // = i(A — 1) is a 
zero of the Wronskian 

W{fo, f+){p) = fo{;p)fi{;p) - a-,p)f+{;p). 

We are mainly interested in eigenvalues A with ReA < 1 and this translates into the condition 
Im^ < since p = — ImA + zRe(A — 1). However, zeros of the Wronskian VF(/o, with 
Im/i < are precisely the resonances or scattering poles of the Schrodinger operator defined by 
Eq. dam). 

4. The linear perturbation theory 

4.1. First— order formulation. We intend to rewrite the wave maps equation (12.50 as a system 
which is first-order in time. From now on we are a bit more careful and specify the full Cauchy 
problem we want to study, including the initial data and the domain. Recall the definition of the 
backward lightcone 

Ct :={(t,r) :tG (0,r),r e [0,T-i]}. 
Precisely, we are interested in the evolution problem 

Mt,r) - Ar{t,r) - ^Mt,r) + ^M^^ = for {t,r) G Ct 
^{0, r) = /(r), Vt(0, r) = g{r) for r G [0, T] 



(4.17) 



with given initial data {f,g)- For the remainder of this section we assume that a solution ^ of the 
Cauchy problem (|4.17p exists and it is sufficiently regular so that all of the following computations 
are justified. In particular, we assume that ip satisfies the regularity condition ^(t,0) = which 
implies /(O) = g{0) = 0. Actually, we want to study the flow about the self-similar solution ■ip'^ 
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and thus, we are interested in "small" perturbations ip of ilP" . Thus, it is convenient to rewrite the 
problem relative to iji^ by setting ijj = ijj^ + Lp and expanding the nonlinearity in Eq. (j4.17p as 

sin(2('(/'^ + p)) = sin(2V''^) + 2 cos(2'(/'^)v? + Nrip) 

where Nt{x) = O(x^) for x — )• is the nonlinear remainder. The Cauchy problem (j4.17p becomes 



(4.18) 



2 I 2 , 2cos(2V'^)-2 
'Ptt - Prr - :;:Pr + ^'P ^ -p-^ P'. 

(/.(O, r) = /(r) - ^'^(O, r), (/.^(O, r) = g{r) - V^RO, r) for r G [0, T] 



NtM. 



in Ct 



where we have suppressed the arguments of if and i/j"^ so as not to overload the equation with 
too many symbols. From the regularity requirement 0) = we obtain the boundary condition 
if{t,0) = for all t G [0, T). Next, motivated by the discussion in Sec. O we introduce the new 
variable (p{t,r) := r'^ip{t,r) which yields 



(4.19) 



'ftt - ^rr + Ih^t + ) + Ar^(r = in Ct 

(^(0,r) =r2[/(r) - V^(0,r)],(^t(0,r) = r2[c,(r) - VHO^ ^)] for r G [0,r] 



and again, we occasionally omit the arguments of and V'^- In order to write this Cauchy problem 
as a first-order system in time, we introduce the variables 



V9i(t,r) : = 



T-t ■ - r 
This choice is motivated by the fact that ipj, j = 1,2, scale like the original field (/? or, in other 
words, they have the same physical dimension. Furthermore, the higher energy (12. Sp is a simple 
expression in terms of ipj. Note also that the field (f can be reconstructed from Lp2 by 



ip{t,r) 



r'ip2{t, r')dr' 



thanks to the boundary condition <p{t,0) = for all t G [0, T). We obtain the system 

dt^2 = ^dr(pi 



(4.20) 



in Ci 



VPi(0,r) = ir2 [g{r)-i^fiO,r)] 
(^2(0, r) = r-'^dy [/(r) - ^'^(0, r)] 



for r G [0, T] 



where J r(p2 is shorthand for r'ip2{t,r')dr' . Now we switch to similarity coordinates 

r = -log(r-t), p=j^^. 

Then we have 



r'(p2{t, r')dr' 



'2t 



(4.21) 



log T. Thus, the system (I4.20p transforms into 




and t = corresponds to r 
dr(l)i = -pdpq 

<9r02 = ^dp(j)l - pdp(j)2 

M- log T, p) = Tp" [g{Tp) - ^IjfiO, Tp)] 

M- log T, p) = Tp [f'{Tp) - i^JiO, Tp)] + 2 [f{Tp) - ^'^(O, Tp)] 

where (/>j (r, /)) := ipj{T — e~'^ ,e~'^ p), j = 1,2 and fo{p) = 2arctan(/9) is the fundamental self-similar 
solution in similarity coordinates. Tracing back the above transformations, we see that {(f)i,(p2) is 
a (sufficiently regular) solution of (|4.2ip if and only if ^ given by 

r'<P2[-logiT-t),^)dr' 
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^(t,r) = V'^(t,r) + 



is a (sufficiently regular) solution of ()4.17p . From now on we restrict ourselves to the linearized 
problem, i.e., we drop the nonlinearity in Eq. (j4.2ip and consider the system 



(4.22) 



5r0i = -pdpcpi + (1)1 + pdp(l)2 - 02 - ^""""p/"^ ^ / P02 1 

9r(/'2 = -pdp4)l - pdp(j)2 J 

</)!(- logr,p) = V [9{Tp) - i'I{0,Tp)] \ fnrn^ \n 11 

M-logT,p) = Tp [f'{Tp) - i>nO,Tp)] + 2 [f{Tp) - V'^(0,rp)] / ^ I"' '\ 



Our aim is to study Eq. ()4.22p by means of semigroup theory, i.e., we rewrite the problem as an 
ordinary differential equation on a suitable Hilbert space. 

4.2. Function spaces and operator formulation. We set 

Hi := {u G C^[0, 1] : u{0) = n'(0) = 0}, H2 := {u G C7^[0, 1] : n(0) = 0} 
and define two sesquilinear forms (-I-)!, (•|-)2 by 

^ dp 



{u\v)i := / u{p)v'{p)^, {u\v)2 := / u {p)v'{p)dp. 
Jo P^ Jo 

Lemma 4.1. The sesquilinear form (•|-)j defines an inner product on Hj, j = 1,2. 

Proof. Obviously, (•|-)2 is well-defined on all of H2 x H2- Furthermore, for u £ H2, {u\u)2 = 
is equivalent to u' = which, in turn, is equivalent to u = thanks to the boundary condition 
u{0) = 0. This proves the assertion for j = 2. For j = 1 we apply Hardy's inequality to see that 



{u\u),= j^\^^^dp< f\u"{p)\^dp<^ 



P^ ^ Jo 

for u £ Hi which is justified by the boundary condition u'(0) = 0. As a consequence, (-I-)! is 
well-defined on all of Hi x Hi. The rest of the proof is identical to the case j = 2 from above. □ 



As usual, the inner product on Hj induces a norm || • \\j := ^ (-lOji J = 1) 2. Moreover, we denote 
by (•!•) the canonically induced inner product on H. := Hi x H2, i.e., 

(u|v) := {Ui\vi)i + {U2\V2)2 



and II • II := yvTO- The vector space Ti equipped with (-I-) is a pre-Hilbert space and we denote 
hy T-L = Hi X H2 its completion. We have the following convenient density result. 

Lemma 4.2. The set C~(0, 1] = {n e C°°[0, 1] : supp(ti) C (0, 1]} is dense in Hj, j = 1,2. 

Proof. Let u E Hi. By definition, this implies ^ G L^(0, 1) where w{p) := p. It is well-known that 
C^(0, 1) is dense in L^(0, 1) (see e.g., [l], p. 31, Theorem 2.19) and thus, for any given e > 0, there 
exists a f G ^^^(0, 1) such that |||y — ';^||l2(o,i) < ^- Set 

"l^ip) ■= [ pv{p')dp'. 



Then u £ C^(0, 1] and ^^-^ = v{p). Consequently, we obtain ||n — n||i = ||^ — u||2,2(o,i) < e 

and, since Hi is dense in Hi by construction, this implies the claim for Hi. The proof for H2 is 
analogous. □ 

For later reference it is useful to note the following embedding and the fact that the boundary 
conditions for the functions "survive" the completion procedure. 

Lemma 4.3. We have the embedding % ^ C[0, 1] x C[0, 1] and u £ H. implies u(0) = 0. 
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Proof. Note first that for u £ H2 we have 



Hp)\ 



u{p)dp 



1 / /■! \ 1/2 



</ \u'{p)\dp<{l \u'{p)\dp\ = ||n'||i2(o,i) = 11^112 



for all yO G [0, 1] and this computation also implies \u{p)\ < \\u\\i for all u G Hi and p G [0, 1]. By 
approximation, these estimates extend to H2 and Hi, respectively, and we obtain 

11112 II ||2 II ||2 iiii2 

II^IIl°°(0,1)xL°°(0,1) = IFiIIl°°(0,1) + II'"2|Iloo(o,1) < l|u|l 

for all u G by taking the supremum over p G [0, 1]. Now let u G 7^ be arbitrary. By construction, 
there exists a sequence (uj) C % such that ||u — Uj|| — t- as j — )■ cxd. The above inequality implies 



u — u 



■illL°°(0,l)xL°°(0,l) 



< llu-uJI ^0 



Lou{p) :-- 



as J —7- 00 and, since the uniform limit of continuous functions is continuous, this shows that 
u G C[0, 1] X C[0, 1]. In particular, uj — )■ u pointwise and, since Uj(0) = for all j, we obtain 
u(0) = 0. □ 

Now we define a linear operator Lq : 'D(Lq) d % ^ % that represents the spatial differential 
operator in Eq. (j4.22p as follows. We set 

P(Zo) := {u G C3[0, 1] : n(0) = n'(0) = u"(0) = 0} x {u G C^[0, 1] : ^(O) = 0} 

and 

-pu'lip) + Ul{p) + pu'^ip) - U2{p) 
-pU'lip) - pu'2{p) 

It should be emphasized that the definition of T){Lq) is a crucial point in the whole construction. 
One has to be very careful with the boundary conditions so as not to "lose" solutions one is actually 
interested in. To see that this is not the case, recall that, in terms of the original perturbation 93, 
the field ipi is given by 

and by regularity, we may assume that ipt{t,r) = 0{r) as r — )• 0+. Since cpi is nothing but ipi 
in similarity coordinates, we may assume 0i(r, /)) = 0{p^) as p ^ 0+. Analogously, one obtains 
4'2{t,p) = 0{p) as /9 — )• 0+. This justifies the chosen boundary conditions. 

Lemma 4.4. For any u G ^{Lq) we have Lqu G Ti. 

Proof. Let u G T>{Lq). The first component [-LquJi of Lqu is obviously twice continuously differen- 
tiable and satisfies the boundary conditions [Lou]i(0) = [Lou]'^(0) = 0. Thus, we have [Lou]i G Hi. 
For the second component we obviously have [I/ou]2 G C'^(0, 1] and by applying de rHopital's rule 
we obtain 

hm [Lou]2(p) = lim ^= lim = 

p^0+ p^0+ p p-!>0+ 

as well as 

P^0+ /3-i>0+ p'^ P^0+ 2 2 

which shows that [-Lou]2 G H2 as claimed. □ 
As a consequence of Lemma 14.41 Lq can be viewed as an unbounded operator 

Lq : v{U) cn^n 
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and Lemma 14.21 implies that Lq is densely defined. We also define another linear operator L' that 
represents the "potential term" in Eq. (j4.22p by T>{L') = % and 



V,{p) p'u2{p')dp' 







where we have set 
(4.23) 



Vi{p) :-- 



L'u(p) := 

V{p)-2 2cos(2/o(/9))-2 
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p- 



(l+p2)2' 



cf. Eq. (13. lip . Obviously, we have L'u G H for any u G T^{L') and thus, L' may be viewed as a 
linear operator 

L' : v{L') cn^n. 

Lemma 4.5. The operator L' : T){L') cH ^ Ti is bounded. As a consequence, it uniquely extends 
to a bounded linear operator L' :% ^ %. 



Proof. By Lemma 14.31 we have 



(4.24) 



\u\ 



L2{0,1) 



u{p)\ dp < \W\\2dp = \\u\\2 



for all u G Ho. Now let u ^ Ti. Then we have 



\t' I|2 

\L u 



^i(p) / pU2ip')dp +Vi{p)pU2{p) 



dp 



M 1 



L°°(0,1) 



1^ 
P 



p'u2{p')dp' 



dp 



(0,1) / \Mp)\ dp 



< 



\pU2{p)\'^dp+ \\U2\\1 < \\U2\ 



<l|uf 

by Hardy's inequality and the estimate ()4.24p . Since Ti is dense mT-L, L' can be extended to all of 
Ti by an elementary result of functional analysis. □ 

In what follows we assume that L' is defined on all Ti. A preliminary operator formulation of 
Eq. (|4.22p is given by 

(A 2^) I f^(r) = (^ + L')1>(T)forr>-logr 

^^•'^^ \ cI,(-logT) = UT 

where $ : [— logT, co) ^ Ti and 

V [giTp)-^pT{0,Tp)] 
Tp [f'{Tp) - ijJiO, Tp)] + 2 [f{Tp) - V^^(0, Tp)] 



ut{p) :-- 



4.3. Well— posedness of the linearized evolution. For simplicity, we consider the free evolution 

f $(t) = Lo$(r) for r > - logT 



(4.26) 



$(-logr) = UT 



first. Note that the formal solution to this problem is "<^(t) = exp(TLo)uj'" , however, it is not 
clear how to define the exponential in this context since the operator Lq is unbounded. We apply 
semigroup theory to assign a precise meaning to "exp(TLo)"- 
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Proposition 4.6. The operator Lq is closable and its closure Lq generates a strongly continuous 
one-parameter semigroup Sq : [0, oo) B{T-L) satisfying 

||5o(t)|| <e-|- 

for all T >0. In particular, the Cauchy problem 

4^^{t) = Lo^{t) for T > -\ogT 
$(-logr) = ur 

for ut € T^{Lq) has a unique solution which is given by 

$(t) = 5o(r + logr)ur 

for allT> -logT. 

Proof. Let u G T>{Lq). In the following calculation, all integrals run from to 1 and for simplicity 
we omit the arguments of the functions. Integrating by parts, we obtain 



Re(Lou|u) = Re I ^ [—pu'i + «i + p«2 — U2]' u[ +Re J 



1 / 

-Ul - pU2 



Uo 



= -Re j ^u'(u[ + ^U2u[ 

— Re ■^u'iu'2 + Re J -u'lu'2 — Re ^ pii2«2 — / \u'2Y 

1 



1 \<{pr 

2 p 



1 



1 



2./ 



U'2{p)»'l{p) 



Re / —U2U1 + Re / —^U2Ui 



1 



Re / -^UiU2 + Rc / -U1U2 — -|«2(1) 
J P J P 2 

1, 



by Cauchy's inequality. Note carefully that the apparently singular boundary terms at p = vanish 
by de I'Hopital's rule thanks to the boundary condition u'((0) = 0. 
Now let f G C~(0, 1] X C~(0, 1] and define 

/(p):=^+/2(.)- 



Then / G C^{0,1]. Set 



Mp) ■= 



P 



^_^2 , fiP'W- 



Obviously, U2 G C°°[0, 1) and 1x2(0) = 0. Expanding / around p = 1 yields 

/(p) = /(I) - /'(1)(1 - p) + /^(l - p)2 + R^p) 

where the remainder R G C°°[0, 1] satisfies the estimate |i?(p)| < C(l — p)^ for all p G [0,1]. 
Performing an analogous expansion for /' and /", one sees that \R^''\p)\ < Cfe(l - pf-^ for 
= 0, 1, 2, 3 and all p G [0, 1]. Thus, for R{p) := R{p')dp', we have the estimate 



d'' R{p) 
dp'' 1 — p 



3-fe 
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for k = 0,1,2,3 and p € [0, 1] by the product rule. This shows that U2 can be written as 



1 + P 



6 ' l-p 

which implies U2 G C'^[0, 1]. Furthermore, we set 

rp 

Mp) ■= p'^Mp) - p [Mp) + Mp')] dp' 



and the above implies ui £ C^[0, 1]. Clearly, we have ui(0) = n'^(O) = u'/(0) = and putting 
everything together we see that u = {ui,U2) £ T>{Lq). By straightforward differentiation one 
verifies that 

(l-Zo)u = f 

and, since f was arbitrary, we conclude that the range of (1 — Lq) is dense in T-L by Lemma 14.21 
The claim now follows from the Lumer-Phillips Theorem, see [15], p. 83, Theorem 3.15. □ 

Remark 4.7. Note that, strictly speaking, we have not solved the free equation (|4.25p but a slightly 
more general problem where Lq is substituted by its closure Lq. This is typical for the approach 
to PDEs with operator methods. 

Remark 4.8. The operator Sq{t) is bounded for any r > 0. Therefore, we can apply it to general 
elements of H and not only to elements in the subspace P(Lo). By doing so, we obtain more general 
"solutions" of the free equation. Although those "solutions" do not actually solve the free equation 
(they solve it only in a weak sense), they are called mild solutions, cf. [15], p. 146, Definition 6.3. 

We immediately obtain a result on the structure of the spectrum ct(Lo) which will be useful later 
on. 

Corollary 4.9. For the spectrum (j{Lq) of the generator Lq we have 

a(Lo) C {A G C : ReA < -i} 
and the resolvent of Lq satisfies the estimate 



ReA + i 



for all X C with ReA > — ^ 



2 • 



Proof. This follows from the growth estimate ||'S'o(t)|| < e a"^ in Proposition 14.61 and |15] . p. 55, 
Theorem 1.10. □ 

By a standard result from semigroup theory we can now conclude the well-posedness of the 
linearized problem. 

Corollary 4.10. The operator Lq + L' generates a strongly continuous one-parameter semigroup 
S : [0, oo) — )■ B{'H) satisfying 

\\S{T)\\<eM + \\L'\\)r 

for all T >0. In particular, the Cauchy problem 

' |:cD(r) = {Lq + L')^{t) forT>- logT 

$(-logT) = UT 

for Ut € DILq) has a unique solution which is given by 

$(r) = 5(r + logr)uT 

for all T > — log T. 
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Proof. This follows from the bounded perturbation theorem, see [15], p. 158, Theorem 1.3. 



□ 



4.4. Spectral analysis of the generator. We analyse the spectrum of the generator Lq + L' in 
more detail which will eventually allow us to sharpen the growth estimate given in Corollary 14.101 
A crucial observation in this respect is the fact that the perturbation L' is compact. 

Lemma 4.11. The operator L' :% ^H. is compact. 

Proof. We define auxihary operators K on C~(0, 1] and A, M on C^{0, 1] x C^{0, 1] by 
{Ku){p) := j\'u{p')dp\ := ( ) , Mu := ( "^^^i 



Then we have (cf. Lemma 

(4.27) Puf = \\Ku,\\l = f ^^^4^^/^=11-2111.(0 1) < 11^^111^(0 1) = < ||uf 

Jo P 

and, by using Hardy's inequality, 

IIMull^ = ||V,„.||f = /' \yMM»)^ VdpWp)\\^ 
Jo P 



< 11^/111^(0,1) / '-^dp+imu,,, ^-^dp 

Jo P Jo P 



) P^ 

< [ \u[{p)\^dp+\\ui\\i 

Jo 



<l|uf 

for all u G C^(0, 1] x C^(0, 1]. Consequently, by Lemma |M1 K, A and M extend to bounded 
operators K : H2 — )• Hi, A : % ^ % and M : % %. By construction, we have L' = MA and 
therefore, by an elementary result of functional analysis, it suffices to prove compactness of A. 

Now recall from the proof of Lemma 14.31 that ||w||j^i(o.i) ^ ll'^^lb for all u G H2 and suppose that 
(uj) C W is a bounded sequence. This implies that (m2j) is a bounded sequence in H2 and it follows 
that («2j) is a bounded sequence in //^(0, 1). By the compact embedding iJ^(0, 1) CC L^(0, 1) 
we conclude that {u2j) has a subsequence which converges in L2(0, 1). Consequently Eq. (OT]) 
implies that (^u-, ) has a convergent subsequence in % which proves the compactness of A. □ 

Lemma 14.111 implies that adding the perturbation L' only affects the point spectrum as will be 
shown next. From now on we write L := Lq + L' . 

Lemma 4.12. If X £ a{L)\a{Lo) then A G ap{L). 

Proof. The proof is actually a standard argument but in order to keep things as self-contained as 
possible, we present it in full detail. Let A G cr(L)\cr(Lo). The well-known identity 

X-L = {1-L'Rl,{X)){X-Lo) 

implies that 1 G a{L' Rlq{X)) since otherwise the operator X — L would be bounded invertible which 
contradicts A G cr{L). However, as a composition of a compact (see Lemma l4.1ip and a bounded 
operator, L'Rl^^^X) is compact and it follows that 1 G ap{L' Rl^^{X)) by the spectral theorem for 
compact operators (see e.g., [19], p. 185, Theorem 6.26). As a consequence, there exists a nonzero 
{ £ n such that L' RL^,{X)i = f. Set u := RLa{X)i. Then u G V{Lo) = V{L), u / 0, and 
(A — L)u = which shows that A G (Tp{L). □ 
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Before we can study the spectrum of L in detail, we need one more technical lemma. The 
problem is that, strictly speaking, we do not know how the operator L acts on general u G T^{L). 
We only know that its action on u G D{Lq) is given by (Lq + L')u. However, the operator Lq (and 
consequently L) is given by the closure of Lq which is a priori an abstract object. The following 
result closes this gap. 

Lemma 4.13. Let u G V{L). Then u G C^{0,1) x C^{0,1) n C[0,1] x C[0,1], u(0) = and 
(A — -L)u = f for f G "H implies 

ui{p) = p'u2{p) + (A - 2) r p'u2{p')dp' - r p'f2{p')dp' 

Jo Jo 

and 

(4.28) -(1 - p^)u"{p) - 2^—^u'{p) + 2\pu'{p) + A(A + l)u{p) + ^n(/>) 

p p^ 

=^+f2{p)+^ rp'f2{p'w 

P^ Jo 

for all p G (0, 1) where u G C^(0, 1) is defined by 

p 



1 fP 

^(p) ■= — / p'u2ip')dp'- 
P Jo 



(4.29) 



Proof Let u G V{L) = V{Lq). Since ueU, Lemma HJ implies u G C[0, 1] x C[0, 1] and u(0) = 0. 
By definition of the closure, there exists a sequence (uj) C 'D{Lq) C C"'^[0, 1] x C^[0, 1] such that 
Uj — )• u and LqUj — )• Lqu in %. This implies (1 — Lo)uj — )■ (1 — Lo)u in ^ and we set fj := (1 — Lo)uj 
and f := (1 — Lo)u. Then i £ T-L and fj- — )■ f . Explicitly, (1 — Lo)uj = ij reads 

pu'ijip) - Pu'2jip) + U2j{p) = fijip) 

-^Kjip) + Pu'2j{p) + U2j{p) = f2j{p) 

for all yO G [0, 1]. Multiplying the second equation in ()4.29p by p^ and adding the result to the first 
equation we obtain 

-p{l - p^Kjip) + (1 + P^)u2j{p) = fl,{p) + p'f2M- 

According to Lemma 14.31 convergence in % implies uniform convergence of the individual com- 
ponents. Therefore, U2j — ?• ^2, fij — )• /i and f2j — ?• /2 in L°°(0, 1) and thus, the above equation 
implies that u^j converges in L°°{a,h) for any [a, C (0,1). Note that we cannot include the 
endpoints due to the singular factors in front of u^j- Consequently, both U2j and Uj^- converge in 
L°°{a, b) and by an elementary result of analysis, this implies that the limiting function U2 belongs 
to C^[a,b] and — )■ Since [a,b] C (0,1) was arbitrary, we obtain U2 G C-'^(0, 1) and the first 
equation in ()4.29p immediately implies ui G C^(0, 1) as claimed. 

As a consequence, the operator L acts as a classical differential operator on u G L){L) and 
(A — L)u = f for an { £ Ti implies 

pu[{p) + (A - l)ni(p) - pu'^ip) + U2{p) + Vi{p) jPp'u2{p')dp' = flip) 
--p<ip) + pu'iip) + ^Mp) = f2{p) 

for all p G (0,1). Thanks to the boundary condition tii(O) = 0, the second equation in (j4.30p 
necessarily implies 

ni(p) = A2(p) + (A-2) r p'u'2{p')dp' - r p'f2{p')dp' 

Jo Jo 
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and inserting this into the first equation in (j4.30p we obtain 

-p{l - p^Kip) + U2ip) + (2A - 1)pWp) + (A - 1)(A - 2) / p'u2{p')dp' 





+ vi{p) r p'u2{p')dp' = hip) + p'hip) + (A - 1) r p'f2{p')dp'- 

Jo Jo 

Defining u G ^^(0, 1) by 

1 fp 

U{p) ■■= -2 / p'u2{p')dp' 

P^ Jo 

we see that u satisfies Eq. K28\i where we use V{p) = p^Vi{p) + 2, cf. Eq. g23|). □ 

Now we are ready to establish the connection between the spectrum of the generator L and the 
mode stabihty analysis from Sec. [3l 

Proposition 4.14. Suppose A G (^p{L) and ReA > — ^. Then there exists a nonzero function 
u{-,X) £ C°°[0,1] that solves Eq. 

Proof. Let A G o'piL) with ReA > — ^ and suppose u G T>{Lq) is an associated eigenfunction. 
According to Lemma 14.131 (A — L)u = implies 

(4.31) ui{p) = p\2{p) + (A - 2) [' p'u2{p)dp' 

and 







-(1 - p')u"{p) - 2^—^u'ip) + 2Xpu'ip) + A(A + l)n(p) + ^u{p) = 
P P^ 



for 



1 

'^ip) '■= — / p'u2ip')dp'- 
P Jo 



Consequently, u G C^(0, 1) satisfies Eq. (|3.12p for p G (0, 1). Note that Eq. ()4.3ip also implies that 
U2 (and therefore u) is nonzero because otherwise, we would have u = 0, a contradiction to the 
assumption that u is an eigenfunction. It remains to show that in fact u G C°°[0, 1]. By standard 
ODE theory it follows immediately that u G C°°(0, 1) since the coefficients of the equation are 
smooth on (0, 1). Furthermore, by Frobenius' method (see e.g., [30], §6.2) we know the asymptotic 
behavior of solutions to Eq. (I3.12p . Around p = 0, the analytic solution behaves like p and the 
nonanalytic one like p~^, see Sec. [3l However, by de I'Hopital's rule we obtain 

lim u(p) = lim '^^ ^ — "^[^ ^ = lim ^^^^ = 
P^o+ p->o+ p^ p^o+ 2 

due to the boundary condition ^2(0) = 0, see Lemma [4.31 and this already implies that u is analytic 
around p = which shows u G C°°[0, 1). Furthermore, around p = 1, the analytic solution behaves 
like 1 and the nonanalytic one like (1 - py~^, at least if A ^ {1,0, -1, -2, . . . } which we shall 
assume for the moment. Now suppose u is not smooth at p = 1. This implies that u can be written 
as 

00 

u{p) = {l-p)'-'Y.a,il-py+h{p) 
j=0 

for suitable constants Oj G C, j G No, and a suitable function h which is analytic around /? = 1. 
Furthermore, ao 7^ and the series has a positive convergence radius (see [30], §6.2). This shows 
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in particular that n"(p) ~ (1 — p) ^ ^ as p — )• 1— . Note that u has to satisfy the integrabihty 
condition 

(4.32) J' (p) dp = J' \u'2{p)fdp<\\u2\\l<oo 

for any e € (0, 1) where p{p) := p. However, we have {-{p^u)')' {p) ~ (1 — p)^^^'^ and this imphes 



-yu)'] (p) 



V 

dpc^ l' |l-pr2-2ReA^^^^ 
'1-e 



Since we assume ReA > This is a contradiction to (14.32P and therefore, we must have u € 
C°°[0, 1] provided that A ^ {0, 1}. If A E {0, 1} and u is not smooth at p = 1, Frobenius' method 
imphes that u is of the form 

oo oo 

u{p) = J2 - py + c(i - p)'-' iog(i -p)Yl - py 

i=o i=o 

for suitable constants aj,bj,c £ C, j £ Nq, where oq / 0, 6o 7^ and the series have positive 
convergence radii (cf. |30) . §6.2). Since we assume that u is not smooth at p = 1, we must have 
c 7^ and this implies u"{p) ~ (1 — p)~^~''* as p — )• 1— which, as before, contradicts (|4.32p . 
Therefore, in any case, we conclude that u £ C°°[0, 1] which finishes the proof. □ 

Now we are ready to obtain a sufficiently accurate description of the spectrum of L. To this end, 
we define the spectral bound for the mode stability problem Eq. (|3.13p as follows. 

Definition 4.15. The spectral bound sq is defined as 

So := supjReA : Eq. (|3.13p has a nontrivial solution v{-, A) G C°°[0, 1]}. 

We say that the fundamental self-similar solution tp'^ is mode stable iff sq < 0. 

Note that Proposition 13.21 implies that sq < ^. 

Lemma 4.16. For the spectrum cr{L) of L we have 

a{L) C {A G C : ReA < max{-i,so}} U {1} 

where sq is the spectral bound, see Definition \4-15\ Furthermore, 1 G <^p{L) and the associated 
(geometric) eigenspace is one-dimensional and spanned by 

f ^ f 2p^ 

sip) ■- 



(l+p2)2 V P(3 + p2) 

Proof. Set M := {A G C : ReA < max{-i,so}} U {1} and let A G a{L). If ReA < -i then 
ReA < max{— ^,so} ^-nd we have A G M. If ReA > — ^ then Corollary 14.91 implies that A ^ a{Lo). 
Thus, according to Lemma [4. 121 A G crp{L). Consequently, Proposition 14. 141 shows that there exists 
a nontrivial u{-,X) G C~[0, 1] that solves Eq. (f3?T2]) . If A = 1 then trivially A G M. Thus, we 
assume A 7^ 1 . Then Lemma 13.11 implies that Eq. (j3.13p has a nontrivial smooth solution and by 
definition of sq we conclude that sq > ReA > — ^- This shows that ReA < sq = max{— ^,so} and 
therefore, A G M. 

We have g G T){Lq) which implies Lg = (Lq + L')^ and by straightforward differentiation one 
obtains (1 — Lq — L')^ = 0. This shows that 1 G crp{L) and g is an associated eigenfunction. Now 
suppose g is another eigenfunction of L with eigenvalue 1. Invoking Lemma 14.131 we see that 

(4.33) ~gi{p) = p'Up) - r p'Up')dp' 

Jo 
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and 

(4.34) - (1 - p\g'\p) - 2^-^~g'{p) + 2\p~g' {p) + A(A + mp) + ^~g{p) = 

p p^ 

where g{p) := -ij Jq p'g2{p')dp' . Furthermore, g G C°°[0, 1]. However, by Frobenius' method we 
know that the only nontrivial solution to Eq. (14.34P in C°°[0, 1] is a multiple of 

9{p) ■= ^2 f P92{p')dp' = ^ 2 
Jo 1 + 

and consequently, there exists a c G C such that 

— / p'Mp') = 9{p) = C9{p) = — / P'92{p') 
P Jo P Jo 

which implies 52 (p) = cg2{p)- Inserting this in Eq. (j4.33|) we obtain gi{p) = cgi{p) and therefore, 
g = eg. This proves that the geometric multiplicity of the eigenvalue 1 G crp(L) is 1. □ 

4.5. Projection on the stable subspace. As already discussed in SecO the unstable eigenvalue 
1 G <yp{L) is a direct consequence of the time translation symmetry of the original problem. In 
this sense, the instability is artificial and we are actually interested in the time evolution with this 
gauge instability removed. By Lemma |4.16| the eigenspace associated to the eigenvalue 1 G crp{L) 
is spanned by a single vector (the gauge mode) which henceforth will be denoted by g G T^{L). We 
intend to construct a projection that removes the gauge instability. However, we cannot simply use 
an orthogonal projection because this is not compatible with the time evolution (the operator L is 
not self-adjoint). Instead, we need to construct a projection that commutes with the semigroup 
S{t). This can be done in a completely abstract way. To this end, let us recall some well known 
facts from spectral theory, see e.g., [19]. At this point we have to warn the reader that our definition 
of the resolvent -Rl(A) differs from [19] by a sign, i.e., we have Rl{^) = — -^)~^ whereas in [19] 
the convention i?_L(A) = (L — A)~^ is used. This has to be taken into account in the following. 

First of all, the resolvent Rl is a holomorphic functionHon p{L). This follows from the expansion 
of Rl{^) in a Neumann series, see [19] . p. 174, Theorem 6.7. According to Lemma 14.16^ the gauge 
eigenvalue 1 is an isolated point in the spectrum of L which implies that Rl has an isolated 
singularity at 1. Therefore, we can find a curve F which lies completely in p{L) and encloses the 
point 1. For instance, we may simply take a circle with radius | and center 1 in the complex plane, 
see Fig. 14.21 Consequently, by a well-known result from complex analysis, Rl can be expanded in 
a Laurent series 

00 

Rl{X)= 

fc=— 00 

around 1 where the "coefficients" (which are actually operators here) are given by the formula 

^k = ^. I (x-ir'-^RLWdx. 

This expression makes perfect sense since the integral is taken over a continuous function (with 
values in a Banach space, though) and it may be defined simply via Riemann sums, just like in the 
scalar case. Consequently, each is a bounded linear operator on %. From this formula and the 



■^Of course, this function has values in the Banach space B{'H) of bounded hnear operators on T-L. However, it 
turns out that this is no obstacle to developing complex analysis just along the same lines as for functions taking 
values in C, see [19) . 
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Figure 4.2. The structure of the spectrum cr{L). According to numerics, the hghter 
shaded region belongs to the resolvent set p{L). However, this has not been proved 
rigorously so far. On the contrary, the darker shaded region contains the spectrum 
of the free generator Lq . 

resolvent equation -Rl(A) — Rhit^) = (a* — ^)RLWRLi^J'), valid for all A, G piL), it follows that 
_i, see [19], p. 38-39. In other words, P := A_i is a projection. Explicitly, we have 

P = ^. f RL{\)d\. 

Furthermore, from the resolvent equation above and this formula, it follows immediately that P 
commutes with Rl{^) for any A G p{L). This implies PL C LP via ^19j, p. 173, Theorem 6.5. We 
define the subspaces M. := PH. and M := {1 — P)'H. Note that both M. and M are closed since 
M = ker(l — P) and N = kerP and P is bounded. This shows that Ai and N equipped with 
(•|-) are Hilbert spaces. However, the most important feature of the projection P is the fact that 
it decomposes the operator L in the following way. Let Lj^ be defined by V{Lj\f) := V{L) CiM 
and Lji/u := Lu. Since P is idempotent, we have for u G 'D{Lj\f) that u = (1 — P)u and thus, 
Lu = L(l — P)u = (1 — P)Lu G M. Hence, Lj\f may be regarded as a linear operator on the Hilbert 
space M. Of course, an analogous statement is true for Lj^ on Ai and from [19], p. 178, Theorem 
6.17, we have (t(L_a/() = {1} and (t{Lj^) = a{L)\{l}. Loosely speaking, Lj^ is "the operator L with 
the gauge instability removed" . Consequently, we call M the stable subspace. 

Lemma 4.17. The operator L^ : 'D{Lj\f) C A/" — )• generates a strongly continuous one-parameter 
semigroup Sj\f : [0, oo) — )• Bi^M) satisfying 

||5.^(t)|| <e(-Hll^'ll)- 

for all T > 0. 

Proof. Note first that Lj\f is densely defined. To see this, let f & M C Ti he arbitrary. Since 
P(L) is dense in T-L, it follows that there exists a sequence (uj) C L){L) such that Uj — t- f in ?^ as 
j — )■ oo. Furthermore, we have PL C LP which in particular implies that PT){L) C ^^{L) and thus, 
(Puj) cMnV{L) = V{Lj^). By the boundedness of P, we obtain ||Puj - Pf|| < ||uj - f|| 0, 
i.e., Puj — )■ Pf = f in ?^ and hence, Puj — )• f in as j — )• oo. Thus, T){Lj^) is dense in N as 
claimed. 

Moreover, let (uj) C T){Lj^f) be such that Uj — )■ u in and Lj\fUj — )• f as j — )■ oo for some 
f G M. This implies uj — )• u, Luj — f as j — t- oo in and by the closedness of L, we obtain 
u G L)(L) and Lu = f . Consequently, we have u G L){L) n A/" = T>{Lj^) and Lj\/u = Lu = f which 
shows that Lj\f is closed. 
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Now recall the estimate from the proof of Proposition 14.61 which readily implies that 
(4.35) Re(Lu|u) = Re(Lou|u) + Re(L'u|u) = Re(Zou|u) + Re(L'u|u) < -^||uf + ||L'u||||u|| 

.(-^,|.',l),|uf 

for all u G T){Lq). Let u € T^{L)- By the definition of the closure, there exists a sequence 
(uj) C T){Lq) such that Uj — )■ u and Luj — t- Lvl as j — t- oo vhH.. Consequently, 

|(Luj|uj) — (Lu|u)| = |(Luj|uj — u) + {Luj — Lu|u)| 

< ll-Z^Ujll ||uj — u|| + \\L\yj — Lu|| ||u|| — t- 

as J — )■ oo which shows that 

1 iiTdAii ii9 / 1 i|2 



(Lu|u) = lim (Luj|uj) < lim h ||L'|| \\vLj\\ = h \\L'\\ ||u 

and thus, the estimate (I4.35P holds in fact for all u G T^{L). In particular, we have 

Re(LAru|u) = Re(Lu|u) < ^-^ + ||uf 

for all u S T){Lj\j-). Finally, for any fixed A > 1 there exists the resolvent Rl{^) (Lemma 14. 16p and 
thus, if f S AA, we can define u := Rl(\){. By definition of the resolvent, u G T^{L) and, since is 
invariant under -Rl(A) (P commutes with Rl{X), see above), we have in fact u = V{L)r\J\f = V[Lj\f) 
and 

(A - LM)n = (A - L)i?i(A)f = f. 

Consequently, the equation (A — Lj^f)\\ = f has a solution for any f G A/" (provided that A > 1) and 
thus, A — -L_/v is surjective. Invoking the Lumer-Philips Theorem (|15). p. 83, Theorem 3.15) yields 
the existence of the semigroup SV with the stated properties. □ 

Next, we relate the semigroup Sj^ to S. 

Corollary 4.18. The stable subspace M is invariant under S{t) and Sj\j-{t) is the restriction of 
S{t) to M , i.e., Sj\f{T) = S{t)\j^ for all r > 0. 

Proof. Fix To > and let f G be arbitrary. According to Lemma [4. 161 the resolvent Rl{^) exists 
for all T G [0,ro] and all n G N provided that n > tq because in this case we have ^ > ^ > 1. 
Note further that the unique solution u G V{Lj^) of (^ — Lj^)vi = (^ — L)u = f is given by Rl{^){ 
(recall that A/" is invariant under Ri{^)) and thus, Rlj^{^) = RLi^)\Af, the restriction of the full 
resolvent to the subspace M. Consequently, the Post-Widder inversion formula for the Laplace 
transform (see p5], p. 223, Corollary 5.5) yields 



n „ /n 



5Ar(r)f= lim -Rlj^I-) f= lim -Rl{-) f = 5(r)f 

n— >-oo IT \T/J ji^cxD Lr Vt/ 



n „ /n 



n 



for all T G [0, To]. Since tq > and f G A/" were arbitrary, this implies the claim. □ 

4.6. Projection on the unstable subspace. Heuristically speaking, the growth behavior of the 
semigroup Sj^ on the stable subspace should be dictated by the spectral bound sq. Before turning 
to this, however, we have to deal with the following subtlety: we do not yet know how many 
dimensions the unstable subspace M. has. Ideally, we would like to have A4 = (g), i.e., the unstable 
subspace should be spanned by the gauge mode. However, this does not follow from the discussion 
above. In fact, since the operator L is not self-adjoint, this is not true in general. A priori, it may 
well be the case that dimA/1 > 1 or even dimA/( = c«. The number dimA/( is called the algebraic 
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(4.36) - (1 - p')u"{p) - 2^—P^u'[p) + 2pu'{p) + 2u{p) + ^u{p) = g2{p) =: g{p) 



multiplicity of the eigenvalue 1 G crp{L) (as opposed to the geometric multiplicity, the dimension of 
ker(l — L), which is 1 by Lemma l4.16p . First, we rule out the pathological case dim = oo. 

Lemma 4.19. The algebraic multiplicity of the gauge eigenvalue 1 G crp{L) is finite. 

Proof. Suppose that dimA^ = oo. According to [19], p. 239, Theorem 5.28, this implies that 
1 belongs to the essential spectrum ae{L) of L 0. However, by [19], p. 244, Theorem 5.35, the 
essential spectrum does not change under a compact perturbation and L = Lq + L' where L' is 
compact (Lemma 14. lip . Consequently, L and Lq have the same essential spectrum but, according 
to Corollary 14.91 1 does not belong to the spectrum of Lq and, since o"e(Lo) C o"(Lo), it cannot 
belong to (Te(-Lo) either, a contradiction. Therefore, the assumption dim = oo must be false. □ 

As a consequence of Lemma 14.191 L_m is in fact a finite-dimensional operator and we are left 
with a linear algebra problem. Since a{Lj^) = {1}, the spectrum of the finite-dimensional linear 
operator 1 — Lj^ consists of only. By linear algebra, this implies that 1 — Lj^ is nilpotent, i.e., 
there exists an m G N such that (1 — = 0. The following result is the key to proving A4 = (g). 

Lemma 4.20. Let c G C\{0}. Then the equation (1 — L)u = eg does not have a solution u G T){L). 

Proof. We argue by contradiction. Without loss of generality we may set c = 1. Now suppose there 
exists an u G T){L) such that 

(1-L)u = g. 

According to Lemma 14.131 this implies 

""^ \pi - ^— „ , , , 

p p- 

for all p G (0, 1) where 

1 

^(P) ■= — / p'u2ip')dp'- 
P Jo 

As before, de I'Hopital's rule implies that m(0) = 0, see the proof of Proposition 14.14] and we also 
have u G C[0, 1] since U2 G C[0, 1] (Lemma 14. 3p . Explicitly, the right-hand side of Eq. (|4.36p is 
given by 

/ N 9iip) . , . Pj^ + P^) 

9{P) = — + 92{P) = T— 2\2 ' 

p"^ (1 + p'^y 

see Lemma 14.161 In particular, we have g{p) > for all p G (0, 1). A direct computation shows 
that hoip) ■■= and 

l + 9p2 + 6p3iogi^ 

'^^'^ = — mr?) — 

are solutions of the homogeneous version of Eq. (|4.36p . Their Wronskian is 

W{ho,hi){p) = — 277- — 27 

and therefore, the variation of constants formula shows that u must be of the form 

fp fp 
u{p) = coho{p) + cihi{p) - ho{p) / p'^hi{p')g{p')dp' + hi{p) / p'^ho{p')g{p')dp' 

J po J pi 

where cq, ci G C and po, pi G [0, 1] are constants. The boundary condition u(0) = implies 

c^ = - f p'^h^{p')g{p')dp\ 
J pi 



^There are several nonequivalent definitions for tlie essential spectrum of a non self-adjoint operator. We use tfie 
one given in [19) . 
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i.e., the singular terms have to cancel in the limit p — )• 0+ and we are left with 



u{p) = coho{p) - ho{p) r p'^hi{p')g{p')dp' + hi{p) r p'\o{p)9{p)dp' 




Similarly, for p — )• 1— we have hi{p) ~ log(l — p) and thus, hi G L}{2i 1) ai^d we must have 



since u G C[0, 1]. This, however, is impossible since the integrand satisfies p'^hQ{p)g{p) > for all 



Now we can prove the desired result, namely that the unstable subspace is indeed spanned by 
the gauge mode. 

Lemma 4.21. For any m G N we have ker(l — L_m)^ = ker(l — L_m). As a consequence, the 
unstable subspace Ai is spanned by the gauge mode, i.e., Ai = (g), and the algebraic multiplicity 
of the gauge eigenvalue 1 G o'p{L) is 1. 

Proof. Note first that Ai C T^{L). To see this, observe that from [19], p. 178, Theorem 6.17 it 
follows that L_M G ^{-M), i.e., is bounded. In our case, this can also be concluded from the 
fact that dim < oo. Now let u G A^. By the density of ^{L) r\M. in A\ (cf. the proof of Lemma 
I4.17p . it follows that there exists a sequence (uj) C 'F>[L) n Ai such that Uj — t- u in as j — )■ oo. 
Since Lj^ is bounded, we conclude that Luj = L_M\ij — )■ L^^u as j — oo and the closedness of L 
implies that u G T^{L). This proves Ai C T>{L) as claimed. 

Next, we prove the assertion ker(l — Lj^)"^ = ker(l — Lj^) for all m G N. Obviously, this is 
true for m = 1 and we proceed by induction. Thus, we show that ker(l — Lj^)"^ = ker(l — L_m) 
implies ker(l — Lj^)"^^^ = ker(l — Lj^) for an arbitrary m G N. Suppose u G ker(l — L_a4)'""^^, 
i.e., (1 — L_m)'""'"^u = 0. Setting v := (1 — L_a4)™~^u we conclude that (1 — L_a4)^v = and thus, 
(1 — L_m)v must be an element of ker(l — Lj^). Furthermore, since Ai C T>{L), v G T>{L) and by 
Lemma 14.161 we obtain 



for a c G C. However, according to Lemma [4. 20 1 this is only possible if c = and therefore, v itself 
belongs to ker(l — L_vi). By the definition of v, this implies (1 — L_vi)™u = 0, i.e., u G ker(l — L_a4)™. 
Invoking the induction hypothesis, we obtain u G ker(l — Lj^) and thus, we have shown that 
ker(l — L_m)"^~^^ C ker(l — Lj^). The converse inclusion ker(l — L_m) C ker(l — L_a/()™'+^ is trivial 
and we arrive at the claim ker(l — L_m)'^ = ker(l — L_m) for all m G N. 

Finally, we show that Ai = (g). To this end, let u G be arbitrary. Since (1 — Lj^/i) is 
nilpotent, there exists an m G N such that (1 — L_M)^vi = 0, i.e., u G ker(l — L_m)^. By the above, 
this implies u G ker(l — Lj^) and invoking Lemma |4.16|. we infer the existence of a c G C such 
that u = eg. Thus, we have shown that Ai C (g). Consequently, Ai is at most one-dimensional. 
However, A4 cannot be zero-dimensional because this would contradict a{Lj^/[) = {1}. This shows 
that dimA^ = 1, A^ = (g) and the algebraic multiplicity (which equals dimA^ by definition) is 
one. □ 

4.7. The time evolution of linear perturbations. In the previous two sections, we have de- 
composed the Hilbert space % into the direct sum % = Ai®M where the closed subspaces Ai and 
are invariant under the time evolution given by the semigroup S{t). The operator Lj^f, which 
acts on the stable subspace AA, has the same spectrum as L except that the unstable gauge eigen- 
value 1 G o-p{L) is missing. Explicitly, we have a{Lj\f) C {A G C : ReA < max{— i, sq}}, cf. Lemma 
14.161 and the unstable subspace Ai is one-dimensional and spanned by the gauge mode g. This is a 
very satisfactory situation which complies with heuristic expectations. Our aim in this section is to 
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(0,1). 



□ 



(1 - L)v = (1 - Lm)"^ = eg 



obtain estimates on the semigroup Sj\f which describes the time evolution on the stable subspace Af. 
In other words, we want to translate the spectral bounds we know for into semigroup bounds for 
the associated time evolution. Since L_\f is not self-adjoint, this is nontrivial. Luckily, there exists 
an abstract machinery which does most of the work for us, the Gearhardt-Priiss-Hwang-Greiner 
Theorem, see e.g., [15], p. 302, Theorem 1.11, or jl8] . This theorem asserts that, if we can find an 
a; G M such that the resolvent Rlj^{X) is uniformly bounded for all A S C with ReA > u, then the 
semigroup Sj\f satisfies the growth estimate [[/^^^(r)!! < e'^'^ for all r > 0. Naively, one might expect 
the estimate H^atIt)!! < e'"^^^-^*"}^ since a{L^) C {A G C : ReA < max{-i,so}}. Based on the 
Gearhardt-Priiss-Hwang-Greiner Theorem, we will prove that, up to an e-loss, this is indeed the 
case. 

First, we obtain a bound for the resolvent i?Lp(A) of the free operator Lq which is particularly 
useful far away from the real axis. In what follows, we denote by [i?L(,(A)f]j, j = 1, 2, the individual 
components of -R^p (A)f. Furthermore, recall the definition of the bounded operator K : H2 — )• Hi 
from the proof of Lemma 14.111 which is given by 

{Ku){p)= r p'u{p')dp'. 



Lemma 4.22. Let e > be fixed but arbitrary. Then the second component of the free resolvent 
satisfies the bound 

Ifll 



\K[RLoim 



2111 < 



|A-2| 

for alii £ v. and all \ £ C with ReA > -| + e, A / 2. 

Proof Let ReA > — ^ + e. According to Corollary HSJ RLf^{X) exists and we have the estimate 

II^Lo(A)||<— ^ 
ReA + ^ 

which shows that 

(4.37) ||[i?Lo(A)f],lb- < ||i?Lo(A)f|| < -fL 

sxeA + 2 

for j = 1, 2 and all f £ H. By setting Vi = in Lemma |4.13[ we see that (A — Lq)u = f implies 

ni(p) = p\2{p) + (A - 2) r p'u2{p')dp' - r p'f2{p')dp' 

Jo Jo 

and, in terms of resolvent components, this reads 

[i?i„(A)f]i(p) = p^[RL,iX)f]2ip) + (A - 2)K[RL,{X)f]2ip) - Kf2{p). 
Consequently, we have 

|A - 2|^||i^[i?z.o(A)f]2||? < ||[i?Lo(A)f]i||? + \p'u'2{p) + 2pu2{pf % + \\Kf2\\l 

Jo P 



< "^"^ + l\^\v!2{p)?dp^ [\u2ipfdp+\\f, 



2II2 



|ReA + i|2 Jo JO 

<||ff + ||[i?Lo(A)f]2||i+||ff 
< l|ff 

r\j W W 

by Lemma 14.31 and the estimate (I4.37p . □ 

Lemma [4. 221 is crucial in order to obtain estimates for the full resolvent Rl{X) far away from the 
real axis. 
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Lemma 4.23. Let e > be fixed but arbitrary. Assume that ReA > — ^ + e and |A| is sufficiently 
large. Then the resolvent Rl{^) exists and it is uniformly bounded as |A| — )• oo. In particular, the 
operator L does not have spectrum jar away from the real axis in the right half -plane {A E C : 
ReA > -\+e}. 

Proof. Let ReA > -\ + e. By Corollary 1131 we see that A G p{Lo) and therefore, we have the 
identity 

\-L = {1-L'Rl,{X)){X-Lo). 

Now note that 



and thus, Lemma 14.221 implies 



\L'RL,[m\ = \\ViK[Rl,{XW 



< 



for all f & Ti. This estimate implies that the Neumann series 

oo 

{i-L'RL,{x)r' = Y.(^'RLoW)'' 

k=0 

converges in the operator norm provided that |A| is sufficiently large. As a consequence, 

Rl{X) = Rlo{X){1 - L'RL,{X)r^ 
exists for large |A| and is uniformly bounded as |A| — )• oo since 

1 



I^Lo(A)|| < 



ReA + i 



by Corollarv 14.91 In particular, there is no spectrum of L far away from the real axis in {A G C : 
ReA>-i+e}. □ 

With these preparations we can prove the crucial estimate for the semigroup on the stable 
subspace which shows that mode stability of the fundamental self-similar solution implies its linear 
stability. This result concludes the linear perturbation theory. 

Theorem 4.24 (Mode stability implies linear stability). Fix £ > 0. There exists a projection 
P G B{'H) onto (g) which commutes with the semigroup S{t) such that 

||S(T)Pf|| = e^l-Pfll 

as well as 

\\S{t){1 - P)f|| < e("^^^{-i^oH=)^||(l - P)f|| 
for all T > and all i €z Ti where sq is the spectral bound, see Definition \4.15\ 

Proof. The existence of P and the fact that P commutes with S{t) have already been discussed 
in Sec. 14.51 Furthermore, by Lemma 14.211 the unstable subspace M = PTi is spanned by the 
gauge mode g. Thus, for any i £ Ti, there exists a constant c(f) G C, depending on f, such that 

II -Pf 1 1 

Pf = c(f)g. Consequently, |c(f)| = -^^ji^ and applying the time evolution yields 

||5(T)Pf|| = |c(f)|||5(T)g|| = enc(f)|||g|| = e^lPfll 

since g is an eigenfunction of the generator L with eigenvalue 1, see Lemma 14.161 

In order to prove the estimate on the stable subspace, note that the statement from Lemma [4.231 
also holds true for the operator Lj\f since Rlj^{X) is the restriction of -Rl(A) to the stable subspace 
A/" = (1 — P)^.. As already mentioned, the spectrum of Lj^ satisfies cr{Lj^) C {A G C : ReA < 
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max{— ^, So}}- This and Lemma 14.231 imply that the resolvent Rlj^{X) is uniformly bounded for 
ReA > max{ — ^, sq} + e, i.e., there exists a constant C > such that 

\\RLAm<c 

for all A G C with ReA > max{— i,so} + Thus, the Gearhardt-Priiss-Hwang-Greiner Theorem 
(|15j. p. 302, Theorem 1.11, see also the recent p!8]) implies the claim. □ 
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